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Abstract. Carbon and nitrogen austenites, modeled by FegN and FegC superstruc- 
tures are studied by full-potential LAPW method. Structure parameters, electronic 
and magnetic properties as well as hyperfine interaction parameters are obtained. Cal- 
culations prove that Fe-C austenite can be successfully modeled by ordered FegC su- 
perstructure. The results show that chemical Fe-C bond in FegC has higher covalent 
part than in FcgN. Detailed analysis of electric field gradient formation for both sys- 
tems is performed. The calculation of electric field gradient allow us to carry out a 
good interpretation of Mossbauer spectra for Fe-C and Fe-N systems. 

PACS numbers: 71.20.Be, 76.80.+y, 71.15.Ap 
1. Introduction 

Face centered cubic (fee) iron-based alloys are widely used for developing of 
stainless austenitic steels especially for using in critical temperature ranges, aggressive 
environment and other severe external conditions. Doping of these steels with light 
interstitial impurities influence mechanics and kinetics of structure phase transitions in 
Fe-based alloys. Nitrogen doping enables to solve the problem of the strengthening of 
stainless steels. 

Carbon and nitrogen doping differently influence structure and mechanical 
properties of austenites. Studying of this effect by means of Mossbauer spectroscopy 
determined substantial differences in spectra of carbon and nitrogen austenitic steels 
(Sozinov et al 1997, 1999). A lot of information, containing in Mossbauer spectra as well 
as high accuracy of measurements lead to the fact that detailed interpretation of these 
experimental results appears to be a very complicated problem and needs application 
of high-accuracy theoretical methods. 

The role of ah initio calculations in the area of solid state research rapidly increased 
last years. Step by step semiempirical methods, containing a lot of parameters, are 
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replaced by calculations from first principles. We observe the extension of application 
field of ab initio quantum mechanical methods of calculations of electronic structure and 
physical properties of solids. Last years these methods are more often used for solving 
not only fundamental, but also applied problems. 

In this paper we study the influence of carbon and nitrogen on atomic and electronic 
structure of fcc-iron. Mossbauer spectroscopy gives the most interesting data about 
impurity distributions, electronic structure and magnetic interactions in solids. The 
study of shifts and splitting of nuclear energy levels gives information about symmetry 
of charge distributions near the nucleus, electronic configurations of atoms and ions 
as well as about peculiarities of atomic structure of solids. We believe that detailed 
interpretation of these experimental data using single up-to-date ab-initio approach is 
an important step in investigation of the influence of light impurities (C,N) in case of 
real Fe(fcc) - based alloys. 

2. Atomic structure 

Up to now there is no full understanding of influence of carbon and nitrogen atoms 
on atomic and electronic structure of fee Fe-C and Fe-N alloys. Based on Mossbauer 
spectroscopy data different authors made different conclusions about atomic structure of 
carbon and nitrogen Fe-based fee alloys. For example, the group of Genin concluded that 
carbon distribution in fee Fe-C alloy is close to the ordered FcgCi-^ structure (Bauer et 
al 1990). On the other hand, the group of Gavriljuk using another method of Mossbauer 
spectrum deconvolution and after performing Monte Carlo computer simulations showed 
that formation of the ordered Fe8Ci_-r structure is improbable (Sozinov et al 1997). 

In this paper we are trying to solve two problems: to study differences in electronic 
structure of carbon and nitrogen austenites, and to perform detailed interpretation 
of Mossbauer experimental data for Fe-N and Fe-C alloys. As an object of our 
investigations we chose model FegC and FegN fcc-type superstructures. The same type 
of C and N ordering in this structures allow us to determine peculiarities of electronic 
structure and hyperfine interaction parameters, connected only with different influence 
of two types of impurity atoms. 

Calculations were performed by applying highly accurate WIEN97 programme 
package (Blaha et al 1999), which employs full-potential LAPW method (Singh 
1994a). Generalized gradient approximation (GGA) according to Perdew-Burke- 
Ernzerhof (Perdew et al 1996) model was used for exchange-correlation potential. The 
value of plane-wave cutoff was Rmt x K„iax = 8.4 which corresponds to about 190 
plane waves per atom in the basis set. Inside atomic spheres the wave function was 
decomposed up to Imax = 12. Charge density and potential was decomposed inside 
atomic spheres using lattice harmonics basis up to L^nax = 6- In the interstitial region 
Fourier expansion was used with 847 coefficients. Calculations were performed for 3000 
K-points in the Brrillouin zone. The radii of atomic spheres were chosen as 1.9 a.u. 
for Fe atoms and 1.6 a.u for N(C) atoms. Standard LAPW basis set was expanded 
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Figure 1. Ordered FegA superstructures (A=C,N). Atoms Fei have one impurity 
atom in the first coordination sphere. Atoms Feg have no impurity in the nearest 
neighbourhood. 

by including local orbitals (Singh 1994b) for Fe3s3p - states and N(C)2s - states. The 
values of all parameters were tested for convergence and ensure accuracy of 0.1 mRy in 
total energy of the system. In order to check the role of local magnetic moments on Fe 
atoms, all calculations were performed using spin-polarized approximation. 

Figure 1 presents FegA supercell (Fm3m, No. 225), which was used for calculations. 
The structure has two symmetry types of Fe atoms: Fei type forms octahedron around 
impurity atom (has one impurity atom in the first coordination sphere); Fcq type has no 
impurity atoms in the first coordination sphere. Before we performed the final electronic 
structure and hyperfine interaction parameters calculations, total energy minimization 
procedure was performed to obtain equilibrium positions of atoms. The lattice constant 
as well as the size of Fei-octahedron were varied to achieve total energy minimum. 
Optimization results were approximated by second-order polynomials using least square 
fit method. Then we analytically found the values of lattice parameters and positions of 
Fe-octahedron atoms. Calculated and experimental values of lattice parameters as well 
as calculated distance between Fe and C(N) atoms are presented in table |1|. Optimization 
results for iron nitride Fe4N are also presented. Good agreement with experiment shows 
that the method can be successfully applied for iron-nitrogen (iron-carbon) compounds. 

3. Electronic structure 

Figure 2 presents total and partial density of states for Fe, C and N atoms in FegA 
compounds. The main difference in electronic structure of two compounds is the energy 
position of impurity 2p states. In FegC carbon 2p states are located right below the 
(i-band of iron. In FegN nitrogen 2p states are located 2 eV below the bottom of Fe 
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Table 1. Structure and magnetic parameters for FegN, FegC and Fe4N compounds. 
Fe-N(C) is the distance between Fe and impurity atom (a.u.) and n is local magnetic 
moment (/Us). 





Lattice 
parameter 


Fe-A 
distance 








fJ-N.C 


/i(aver./ 
atom) 


FesN (theor.) 


13.946 


3.54 


2.78 


2.09 




-0.11 


2.25 


(expt) 


13.849* 














FesC (thcor.) 


13.981 


3.58 


2.71 


2.12 




-0.16 


2.25 


(expt) 


13.839^ 














Fe4N (theor.) 


7.164 


3.582 


2.90 




2.31 


0.04 


2.50 


(expt) 


7.171'' 


3.585^ 


S.OO'^ 




2.00^^ 




2.21^^ 



Cheng et al (1990). 
^ Jacobs et al (1995). 
" Frazer (1958). 



d-band. This difference in energy positions of C and N 2p states leads to essential 
difference in the character of chemical bonding in these two compounds. In FegC we 
see stronger p — d hybridization compare to FcgN compound. The width of N 2p-band 
is twice shorter than corresponding 2|)-band of carbon. This probably leads to much 
higher Fe-C interaction compare to Fe-N interaction. 

The difference in energy distribution of valent electrons leads to considerable 
difference in electronic density distribution in FegN and FegC. Figure 3 presents 
electronic density distribution in (100) plane for two energy intervals: 2p states of 
impurity atoms and 3d states of iron in FegN and FcgC compounds. Nitrogen 2p 
electrons are much more localized (figure 3a) compare to 2p electrons of carbon (figure 
3c). Charge distribution around Fei atom in FegC is more symmetric (figure 3d) 
compare to corresponding charge distribution in FegN (figure 3b). It should be noted 
that impurity atoms cause asymmetric charge distribution around Fei atom. This also 
leads to asymmetric charge distribution of d electrons around Fei atom. These effects 
greatly influence electric field gradient (EFG) at Fei nuclei, which can be measured as 
quadrupole splitting in Mossbauer spectra of these two compounds. 

4. Hyperfine interactions in FegN and FegC 

For our structures ^''Fe nucleus quadrupole splitting is given by the following expression: 
A = ^eQV,,, (1) 

where Q - quadrupole moment of the nucleus and V^z - principal component of electric 
field gradient (EFG) tensor. 

In our calculations we used Q^^Fe =0.166, which was determined by Dufek, Blaha 
and Schwarz (1995) by comparing experimental quadrupole splitting and calculated 
EFG values for a large number of different Fe compounds. EFG is calculated on ab- 
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Figure 2. Density of states (DOS) for FesA (A=C,N) compounds 

initio basis directly from electronic density distribution using method, developed by 
Blaha, Schwarz and Herzig (1985). 

Table |] presents our theoretical results for quadrupole splitting as well as 
experimental data for carbon and nitrogen austenites, obtained by several groups. We 
see that reasonable agreement with experiment is obtained for Fe-C alloy. This leads to 
conclusion that local environment of Fei atom in real alloy is close to Fei environment 
in our model FegC superstructure. This is also supported by the fact that according 
to experimental data (Bauer et al 1988, Gavriljuk and Nadutov 1983, Oda et al 1994) 
Fe-C austenite has only Feo and Fei atoms, which corresponds to Fe configurations in 
our model superstructure. On the other side, in Fe-N austenite Mossbauer spectrum 
is formed by contributions from Fcq, Fei and Fe2_i8o° (dumb-bell configuration) atoms 
(Oda et al 1990, Foct et al 1987, Gavriljuk et al 1990). Both absence of Fe2-i8o° atoms 
in our structure and bad agreement with experiment shows that Fe-N austenite can not 
be modeled by FesN superstructure. 

Our calculation shows that contribution to Vzz (the principal component of EFG 
tensor) from the regions outside the spheres (lattice EFG) is 5%, 3% and 10% for FegN, 
FcgC and Fe4N respectively. So, for understanding of the origin of the EFG we focus 
on its main component, the valence EFG, which originates from the space inside the 
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(a) FeBN N-2p charge distribution in (1DD) plane, e/A"3 

(c) Fe8C C-2p charge distribution in (100) plane, e/A**3 




(b) FeBN Fe-3d charge distribution in (IDD) plane, e'A"3 

(d) FeBC Fe-3d charge distribution in (1 00) plane, eWS 




Figure 3. FegN and FegC charge distribution in (100) plane (e/A^) for different 
energy regions. 

atomic spheres. The ingredients for the calculation of the valence EFG are the density 
coefficients P2M, which originate from two radial wave functions with / and /': 

PLAiir) = E E E Rim{r)R,n.'{r)GZr\ (2) 

E<Ef l,m l',m' 

where Rim{r) - LAPW radial wave functions, Gf^;™™"' - Gaunt integrals. 

We performed comparative analysis of EFG formation for FegN and FegC structures. 
Table ^ and table § show that total EFG and quadrupole splitting is almost 3 times 
larger in FegC structure compare to FegN. By performing the analysis of mechanisms of 
EFG formation, we tried to find out the reasons of this considerable difference. First of 
all it is obvious that contribution from Fe 3s3p energy interval remains constant while 
going from nitrogen to carbon compound. This means that redistribution of electronic 
density in this region is not considerable and can not infiuence total EFG increasing. 
We observe increase of negative contribution from C 2s interval, but it is also not 
considerable. Moreover, the sum of contributions from impurity 2s2p energy interval 
remains almost constant while going from nitrogen to carbon. This means that the main 
reason of EFG increasing is redistribution of charge density in Fe band (and partially 
in impurity 2p band). 

In FegN 2p3(i-band EFG contribution is positive, which decreases total EFG. In 
FegC this contribution is negative and increases total EFG. In FegC charge distributions 
in 2p and ?)d bands is more symmetric than in FegN. While going from N to C the relative 
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Table 2. Calculated and experimental values of quadrupole splitting A(mm/s) in 



Mossbauer spectra for Fe-N, Fe-C fee alloys. 


Fe4N FesN Feio.iN FenN 


FegC 


Fen.aC Fen.sC 


A (theor.) 0.50 0.18 


0.53 




(expt) 0.50^ 0.39^ 0.25'= 




0.64'^ 0.672 


0.39^^ 




0.63^ 



Foot (1974). 
^ Oda et al (1990). 
" Foot et al (1987). 

Gavriljuk et al (1990). 
" Genin and Flinn (1968). 
f DeCristofaro and Kaplow (1977). 
e Oda et al (1994). 



Table 3. EFG analysis for Fei atom in FcsN and FcsC superstructures. Contributions 
to valent EFG (10^^ V/m^) from different energy bands are presented. 



FegN 


yp-p 

^ zz 


^rd—d 

zz 


\rval 

zz 


FegC 


yp-p 

zz 


'ird—d 

zz 


zz 


Fe3s3p 


+2.448 


0.000 


+2.417 


Fe3s3p 


+2.455 


0.000 


+2.434 


N2s 


-4.369 


-0.174 


-4.554 


C2s 


-4.822 


-0.354 


-5.194 


N2p 


-4.552 


-2.467 


-7.067 


G2p 


-2.479 


-3.177 


-5.695 


Fe3rf 


+0.686 


+7.310 


+8.074 


Fe3rf 


-1.200 


+6.309 


+5.166 


Valent 








Valent 








hand 


-5.787 


+4.669 


-1.130 


hand 


-6.046 


+2.778 


-3.289 



increasing of symmetry of 3(i charge distribution is much larger than relative increasing 
of symmetry of 2p charge distribution, which leads to smaller EFG compensation and, 
as a result, to negative contribution from 2p3d-band. 

Let us now consider partial EFG contributions {p—p, d — d) inside Fci sphere (other 
symmetry contributions like s — d, etc. are negligibly small and can not determine EFG 
formation). While going from N to C we observe relatively small increase of the absolute 
value of p — p contribution and quite large decrease oi d — d contribution. As far as 
d — d contribution is mostly formed by own Fe3d electrons, we can make a conclusion 
that while going from N to C the increasing of space symmetry of Fe 3d electrons causes 
stronger influence of the asymmetric distribution of 2s and 2p electrons of the impurity, 
located in Fci sphere. 

Finally we tried to analyze the reasons of symmetry increase of 3(i-band charge 
distribution in FcgC Let us consider partial EFG contributions only for 3(i-band. In 
FegN p — p and d — d contributions are both positive and decrease total (negative) EFG. 
In FcgC situation is different: p—p contribution becomes negative and its absolute value 
becomes twice larger. This decreases positive d — d contribution and causes the increase 
of total EFG. 
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5. Summary 

The results of our full-potential band structure calculations of FegN and FesC structures 
lead to the following conclusions. We observe that in FegC chemical Fe-C bonding is 
more covalent compare to Fe-N bonding in FegN. Fe-C austenite can be successfully 
modeled by ordered FesC superstructure, containing iron atoms in Feo and Fei 
configurations. Calculated quadrupole splitting is in good agreement with experimental 
values. 

We performed detailed analysis of EFG formation on Fe nuclei in this structure. 
The main EFG contribution in made by carbon 2s and 2p states. The contributions 
of these two states appear to be almost equal. Carbon states contributions are greatly 
canceled by asymmetry of space distribution of FeSci-electrons. 

We also calculated structure, magnetic and hyperfine interaction parameters for 
ordered iron nitride Fe4N structure. The results of our calculations are in excellent 
agreement with experimental data. 

We come to conclusion that Fe-N austenite can not be modeled by FegN structure. 
The results of our calculations show that nitrogen distribution in real Fe-N alloy is 
different from nitrogen distribution in FegN structure. 
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